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EARLY  TIME  COUPLING  STUDIES  USING  A  ID  HYBRID  CODE 
I .  INTRODUCTION 

The  understanding  of  atmospheric  effects  generated  by  a  HANE  is  of 
vital  importance  to  DNA.  It  is  well  known  that  HANE's  cause  severe 

3 

disturbances  in  the  atmosphere  and  ionosphere  which  impact  C  I  systems.  In 
particular,  long  lasting,  large  scale  ionization  irregularities  degrade 
radar  and  communications  signals.  To  develop  predictive  capability  for 
such  HANE  effects,  DNA  sponsored  a  major  research  effort  at  NRL  leading  to 
theoretical  and  computational  models  of  HANE.  Significant  progress  was 

made  in  this  area  during  the  past  decade  leading  to  the  development  of 
sophisticated  multifluid  codes  which  incorporate  in  a  self  consistent 
manner  anomalous  transport  due  to  instabilities.  Recent  advances  in  both 
plasma  theory  and  numerical  techniques  allow  us  today  to  produce 
substantially  more  refined  and  accurate  models  of  the  early  time  HANE 
phenomenology.  These  advances  indicate  the  need  for  a  better  incorporation 
of  kinetic  effects  due  to  strong  deviations  of  the  particle  distributions 
from  Maxwellian.  The  new  codes  allow  us  to  compute  the  laminar  and 

turbulent  dynamic  ion  response  over  long  time  and  space  scales  and  for  a 

mixture  of  ionized  species.  In  this  way  laminar,  electrostatic,  Larmor  and 
turbulent  coupling  processes  can  be  followed  dynamically  and  their  relative 
strength  and  dynamic  interplay  can  be  assessed.  This  paper  is  the  first  to 
describe  this  series  of  investigations,  focuses  on  a  particular  process, 
early  time  debris-air  coupling,  and  stresses  quantitative  understanding  of 
the  underlying  physics. 

Early  time  coupling  had  been  a  controversial  subject  in  the  DNA  HANE 
community  in  the  1970'3.  Two  "opposing"  views  were  maintained:  NRL,  LASL 
and  ARA  advocated  short-range  coupling  associated  with  plasma  turbulence, 
while  MRC  argued  for  "Larmor  coupling"  associated  with  the  ion  gyration  in 
the  ambient  magnetic  field.  To  a  certain  extent,  the  issue  was  not  whether 
or  not  these  coupling  mechanisms  worked,  but  which  mechanism  was 
dominant.  The  turbulent  coupling  process  developed  at  NRL  was  studied  by 
Incorporating  theoretical  estimates  of  anomalous  transport  coefficients  in 
a  ID  multi-fluid  code.  Although  this  procedure  yielded  results  consistent 
Manuscript  approved  February  14, 1985. 
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with  HANE  data  and  provided  strong  support  for  turbulent  coupling.  It  was 
not  a  self-consistent  technique.  Moreover,  It  was  not  capable  of  handling 
kinetic  effects  such  as  caused  by  reflected  ions;  such  reflected  Ions  are 
known  to  occur  In  high  Mach  number  shocks.  The  present  study,  using  a  ID 
hybrid  code,  overcomes  these  deficiencies.  Namely,  It  allows  Instabilities 
to  be  excited  which  affect  the  plasma  evolution  in  a  self-consistent 
manner,  i.e.,  no  model  anomalous  transport  coefficients  are  used,  and 
reflected  ions  are  permitted.  Furthermore,  the  ID  hybrid  code  inherently 
includes  "Larraor  coupling"  so  that  a  direct  comparison  of  the  relevant 
importance  and  interplay  of  the  coupling  mechanisms  can  be  made. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
present  an  overview  of  the  basic  physical  processes  occurring  as  debris 
streams  through  the  air.  In  Section  III  we  present  results  of  two 
simulations  which  demonstrate  the  various  aspects  of  coupling.  Finally,  in 
the  last  section  we  summarize  our  results  and  discuss  their  implications  to 
HANE  phenomenology. 

II.  OVERVIEW  OF  PHYSICAL  PROCESSES 

The  evolution  of  the  debris-air  plasma  is  characterized  by  three 
temporal  stages:  (1)  magnetic  field  compression;  (2)  piston  formation;  and 
(3)  shock  formation  and  evolution.  The  bulk  of  this  paper  describes  in 
detail  piston  dynamics;  for  completeness  we  briefly  discuss  the  other  two 
stages. 

In  the  initial  stage,  the  debris  streams  outward  from  the  burst  point 
and  picks  up  air  electrons  which  merge  with  the  debris  electrons  to  form  a 
single  electron  distribution.  The  compression  of  the  air  electrons  leads 
to  a  compression  of  the  magnetic  field  since  the  air  electrons  are  "tied" 
to  the  field  lines  due  to  the  Low  value  of  resistivity.  This  is  the  first 
stage  of  the  debris-air  interaction  and  has  been  well-documented 
theoretically  (Longraire,  1963;  Sloan,  1970;  Lampe  and  Hernandez,  1972), 
numerically  (Clark  et  al.,  1974),  and  experimentally  (Rlpin  et  al.,  1984). 

The  next  stage  involves  the  dynamics  of  piston  formation  and  will  be 
discussed  in  detail  in  the  following  section.  Briefly,  this  stage  is 
dominated  by  the  interaction  of  the  two  ion  streams,  i.e.,  the  debris  ions 
and  the  air  ions.  It  is  this  phase  where  various  coupling  mechanisms 
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become  apparent.  The  debris  ions  slow  down,  while  the  air  ions  are  picked 
up;  there  is  momentum  exchange  between  the  ion  species.  There  is  also  a 
conversion  of  the  directed  kinetic  energy  of  the  debris  to  thermal  energy 
of  both  ions  and  electrons,  i.e.,  heating  and  acceleration  occurs.  And 
finally,  an  electrostatic  potential  forms  at  the  leading  edge  of  the 
magnetic  field  compression  which  causes  some  debris  ions  to  be  accelerated 
ahead  of  the  magnetic  piston  and  some  air  ions  to  be  reflected.  This  final 
point,  although  well-known,  can  be  very  important  to  the  formation  of 
debris  patches  in  the  conjugate  regions  since  it  results  in  very  high 
velocity  ions. 

Lastly,  the  final  stage  Involves  the  formation  of  a  shock  wave  which 
"runs"  ahead  of  the  piston.  We  do  not  describe  this  phase  in  detail  in 
this  paper,  but  do  show  that  a  minimum  time  is  required  for  high  Mach 
number  shock  formation.  This  minimum  time  is  basically  the  time  needed  for 
the  shock  to  "run"  ahead  of  the  accelerated  debris  or  reflected  air  ions. 


III.  SIMULATIONS  AND  RESULTS 


Prior  to  discussing  the  details  of  the  simulations,  we  present  a 
simple  model  equation  that  elucidates  the  various  coupling  mechanisms  to  be 
studied.  Moreover,  it  highlights  the  distinction  between  previous  models 
(e.g.,  KLYSMA)  and  the  present  work.  We  consider  the  following  momentum 
equation  in  the  x  direction  for  an  air  ion  (i.e.,  radial  direction 

A 

perpendicular  to  the  ambient  magnetic  field  ^ 


eE. 


m 


n  V  -  v*rv  -  V,) 
a  ay  ^  ax  d-^ 


(1) 


where  the  subscript  a  refers  to  air,  D  =  Z  eB/ra  c,  v*  is  an  anomalous  ion- 

A  £13  3 

ion  collision  frequency,  =  ^d®x  streaming  debris  velocity,  e  is 

the  charge,  m  is  the  mass,  and  Z  is  the  charge  state.  The  first  term  on 
the  RHS  of  (1)  arises  from  a  laminar  electric  field  usually  associated  with 
the  leading  edge  of  the  magnetic  compression.  It  acts  to  accelerate  debris 
tons  and  to  reflect  air  ions.  The  second  term  in  (1)  is  the  magnetic  force 
which  is  associated  with  Larraor  coupling.  The  final  term  in  (1) 
corresponds  to  turbulent  "pick  up"  of  the  air  and  arises  because  of  plasma 
Instabilities.  For  plasma  turbulence  such  that  v*  >  it  is  clear  that 
turbulent  coupling  can  dominate  over  Larmor  coupling. 
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Notice  that  the  force  n  V  associated  with  Larmor  coupling  is 
proportional  to  the  value  of  which,  before  any  form  of  coupling  occurs, 


is  V, 


V,.„  =  0.  Therefore  the  laminar  force  eE  /m  and  the  turbulent 
<iy  X  a 


force  -  Vjl  ~  v*Vj  will  dominate  initially  even  for  v*  «  (i.e., 

as  long  as  v*  >  ^gCV^y/Vj)).  The  time  evolution  of  will  be  given  by 

dV  eE  eV.B 

T?-  •  “  V,  (2) 

dt  m  me  ad 

a  a 

In  deriving  eq.  (la)  we  assumed  an  infinitely  conducting  (p  -*•  0),  low  S 
( 8  <<  1)  situation  so  that  E^  =  V^B/c.  From  eqs.  (1)  and  (2  )  it  is 
obvious  that  the  time  scale  for  significant  Larmor  coupling  is  R^t  >  1, 
i.e.,  the  value  of  should  approach  V^.  At  earlier  times  the  laminar 

and  turbulent  forces  will  dominate.  The  turbulent  forces  will  also 
dominate  even  at  t  »  as  long  v*  >  l^ax  ~  '^dl'^'^d*  Another  important 
point  follows  from  (1)  and  (2)  if  we  neglect  the  laminar  and  turbulent 
forces.  If  L  is  the  size  of  magnetic  field  compression  then  as  long  as 
<  Vj  the  maximum  value  of  will  be  V^y  =  mentioned  previously 

the  Larmor  momentum  coupling  force  in  the  x  direction  becomes  significant 
if  V^y  =»  n^L.  Therefore  for  effective  Larmor  coupling  L  >  Namely, 

an  extremely  broad  magnetic  field  profile  will  be  required.  We  will 
further  comment  on  this  point  later  on. 

In  previous  work,  theoretical  estimates  of  v*  were  incorporated  in  a 
multi-fluid  code  (KLYSMA)  which  solved  momentum  equations  similar  to  (1) 
and  (2)  (Clark  et  al.,  1974).  However,  the  turbulent  Interaction  Involves 
ions  interacting  with  fluctuating  electric  fields,  i.e.,  v*  =  v*(6E).  The 
present  simulations  Incorporate  such  effects  self-consistently.  The  air 
momentum  equation  solved,  analogous  to  (1),  is 


-  ^  fE 
m  X 
a 


f  n  V 
a  ay 


where  5E^  is  self  consistently  generated  because  of  plasma  instabilities 
(e.g.,  ion-ion  streaming  instability  (Papadopoulos  et  al.,  1971)  and  is 
allowed  to  act  directly  on  the  ion  motion,  rather  than  through  a  model 
equation  such  as  (1).  Given  this  background,  we  now  present  the  simulation 
results. 


The  simulations  were  performed  using  a  one  dimensional  quasi-neutral 
hybrid  code  (Chodura,  1975;  Wlnske  and  Leroy,  1984)  whose  description  Is 
given  in  Appendix  I.  The  code  Incorporates  ion  kinetic  effects  by 
following  the  trajectories  of  "superparticles",  numerically  representing 
many  actual  ions,  in  3  dimensions  in  velocity  space  and  one  dimension  in 
configuration  space.  The  electrons  are  approximated  as  a  massless  fluid 
described  by  momentum  and  energy  transport  equations;  the  electron  pressure 
is  assumed  to  be  isotropic.  These  equations,  with  Maxwell's  equations,  are 
solved  on  a  one  dimensional  spatial  grid  using  moments  calculated  from  the 
ion  distributions.  Anomalous  resistivity,  resulting  in  magnetic  field 
diffusion  and  electron  ohmic  heating,  is  included  as  a  parameter  in  the 
equations  with  a  value  consistent  with  that  expected  from  current  driven 
cross-field  instabilities. 

The  boundary  conditions  in  the  code  correspond  to  the  conditions  we 
expect  at  early  times.  A  dense  ion  beam  is  injected  continuously  into  the 
simulation  box  which  is  initially  filled  with  a  tenuous  stationary 
plasma.  In  our  present  results  both  ion  populations  are  H'*’,  but  other  ion 
species  can  be  easily  considered  and  will  be  presented  elsewhere.  We  set 
up  a  finite  width  magnetic  field  compression  within  the  leading  edge  of  the 
ion  beam.  This  represents  the  field  compression  expected  from  the  cross- 
field  beam  motion  which,  because  of  finite  resistivity,  will  diffuse  into 
the  beam.  The  width  of  the  compression  has  been  varied  to  study  various 
physical  phenomena.  The  magnetic  field  is  held  fixed  at  the  edges  of  the 
simulation  box. 

We  present  below  the  results  of  two  simulations  selected  to  illustrate 
the  underlying  physical  coupling  processes.  The  initial  state  is  typical 
of  the  piston  formation  stage.  It  is  shown  in  Fig.  la,  for  the  case  of 
=  Vjj/Va  =  4,  where  is  defined  on  the  basis  of  the  ambient  (l.e., 

upstream)  parameters.  A  magnetic  pulse  consistent  with  B/n  =  constant  is 
introduced  at  the  left  hand  boundary.  In  this  paper  we  do  not  study  the 
formation  of  the  pulse  but  assume  its  nature  in  accordance  with  previous 
work  mentioned  in  Section  II,  and  its  width  L  to  illustrate  the  important 
physics.  In  the  figures  the  velocity  is  in  terms  of  the  upstream  value  of 
and  of  time  in  terms  of  the  upstream  value  of  9  The  units  of  length 
are  in  terras  of  the  ion  cyclotron  value  computed  for  air  with  velocity 


but  with  the  value  of  the  downstream  (l.e.,  compressed)  magnetic  field. 
This  representation  has  been  selected  for  clarity  of  the  underlying 
physics.  We  consider  (1)  a  "broad"  Lg  =  (R^  in  upstream  values)  pulse 

with  1  4  and  n^j/n^  =  16  and  (2)  a  "narrow"  Lg  =  pulse  with  =  8 

and  n^/Ug  *  64.  Here  V^,  n^j  are  the  debris  velocity  and  density,  the 
air  density  and  is  defined  on  the  basis  of  the  upstream  conditions. 

A.  Broad  Pulse  (Lg  *• 

This  simulation  is  initialized  with  a  magnetic  pulse  width 
=  4  and  n^/n^  =  16.  We  show  a  series  of  results  for  this  run  in  Fig. 
1.  Figure  1  displays  four  Important  quantities  versus  distance.  The 
bottom  curve  represents  the  magnetic  field  magnitude,  the  middle  curve  (or 
"dots")  represent  particle  velocities  of  the  debris  and  air  ions  (the 
debris  ions  have  ~  4V^  while  the  air  ions  have  ~  0),  and  the  top 
curve  represents  the  electric  field  in  the  x  direction.  Figure  la  is  at  t 
=  0.0059  where  is  based  on  the  upstream  value  of  B.  Note  the  broad 
pulse  in  B  with  sharp  gradients  at  the  leading  and  trailing  edges. 
Associated  with  these  gradients  are  strong  laminar  electric  fields  at  these 
edges  as  shown  in  the  top  portion  of  Fig.  la.  These  fields  strongly  affect 
the  ions.  Note  in  the  middle  panel  that  some  debris  ions  are  being 
accelerated  ahead  of  the  B  pulse,  and  air  ions  are  being  picked  up  in  a 
reflection  process. 

Figure  lb  depicts  the  system  at  t  =  0.0386  and  several  interesting 
features  are  developing.  First,  the  B  field  pulse  has  broadened 

considerably.  The  leading  edge  of  the  pulse  has  been  "dragged"  out  to 
X  ~  4r^^,  and  the  gradients  are  not  as  sharp.  Second,  the  laminar  electric 
field  has  fallen  in  intensity  and  is  accompanied  by  a  considerable  amount 
of  turbulence  in  the  region  0  <  r/r  <  2.4  caused  by  of  the  ion-ion 
streaming.  This  turbulence  causes  some  momentum  exchange  between  debris 
and  air  ions,  and  therraallzation  of  the  particles  in  this  region.  Finally, 
also  at  this  time,  a  fraction  of  debris  ions  have  been  accelerated  to  a 
velocity  greater  than  twice  the  Initial  velocity. 

Figure  Ic  shows  the  system  at  t  =  0.0874  The  magnetic  pulse 

continues  to  broaden  and  now  extends  to  X  ~  br^i.  A  laminar  electric  field 
still  exists  at  x  ~  3.6  However,  the  turbulence  in  the  region  0  < 


x/r^^  <  3.6  has  subsided  since  the  relative  velocity  between  debris  and  air 
ions  is  small.  The  ion  motion  in  this  region  is  dominated  by  the  magnetic 
field  and,  in  effect,  is  where  "Larmor  coupling"  is  taking  place.  Perhaps 
the  most  Interesting  region  is  3.6  <  x/r^j^  <  6.0.  It  is  clear  that  debris 
acceleration  and  air  reflection  are  continuing.  Furthermore,  turbulence  is 
developing  at  the  debris-air  Interface  (x  ~  4.0 

Finally,  Fig.  Id  shows  the  system  at  t  =  0.  1542  0^^.  The  magnetic 
pulse  is  quite  broad  now,  extending  out  to  x  ~  8.0  Very  intense 
electrostatic  turbulence  is  occurring  at  the  debris-air  mix  in  the  region 
4.2  <  x/ r^^  <  8.2  which  is  producing  coupling  between  the  ion  species  and 
thermalizing  the  particles.  Debris  particles  have  been  accelerated  up  to 

V  ~  In  the  region  2.0  <  x/r^^  <  4.0  the  air  ions  have  a  velocity 

V  ~  5.0  V^  while  the  debris  ions  have  a  velocity  V  ~  3.0  V^.  In  this 
region  the  momentum  coupling  is  through  the  magnetic  field  (Larmor 
coupling)  and  the  debris  and  air  ions  perform  gyrations  about  the  magnetic 
field.  Very  little  turbulence  is  evident  in  the  electric  field,  confirming 
the  fact  that  v*  <  0^. 

Thus,  in  this  simulation  it  is  found  that  three  important  processes 
are  affecting  coupling,  energization,  and  thermallzation;  (1)  laminar 
electric  fields;  (2)  turbulent  electric  fields;  and  (3)  magnetic  field 
pick-up.  The  laminar  fields  are  most  prominent  early  in  the  run  (Fig.  la) 
and  cause  debris  acceleration  and  air  reflection.  This  is  crucial  to  early 
time  HANE  processes  since  debris  ions  can  be  accelerated  up  to  ~  3  times 
their  initial  expansion  velocity,  and  air  ions  are  "picked-up".  Turbulent 
aiecttlc  fields  cause  coupling  and  thermaiization  of  the  ion  species.  In 
this  run  they  are  most  effective  at  the  debris-air  Interface,  and  ahead  of 
the  interface  (Figs.  Ic  and  Id).  Finally,  "Larmor  coupling"  is  observed 
within  a  section  of  the  broad  magnetic  pulse  in  Figs.  Ic  and  Id,  i.e.,  the 
debris  and  air  ions  have  rotated  about  the  B  field  so  that  the  air  Ions 
have  a  larger  velocity  than  that  of  the  debris  (V^^^  >  V^^) .  This  process 
supplements  any  incomplete  part  of  the  turbulent  coupling  and  operates  on  a 
long  time  scale. 
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B.  Narrow  pulse  [L-  ==  O.l  r^  .  j 

D  u  1 

This  simulation  is  initialized  with  a  magnetic  pulse  width  =0.1 
and  =  8  and  =  64.  The  results  are  shown  in  Fig.  2.  In 

Fig.  2a  we  show  the  plasma  system  at  t  =  0.0206  Already  by  this  time 

a  substantial  number  of  air  ions  are  being  reflected  from  the  laminar 
electric  field  at  x  ~  1.5  Furthermore,  debris  ions  have  been 

accelerated  up  to  V  ~  2.5  and  extend  out  to  x  ~  this  region 

there  is  a  considerable  amount  of  electrostatic  turbulence  as  seen  in  the 
electric  field  plot.  Figure  2b  is  at  t  =  0.0411  The  magnitude  of  the 

turbulence  electric  fields  have  increased  from  the  previous  time  shown  and 
are  quite  intense  at  the  debris-air  interface  (x  ~  3.0  The  air  ions 

are  being  "picked  up"  because  of  their  turbulence  and  thermalized.  The 

magnetic  field  has  been  dragged  out  to  x  ~  7.0  r^^  by  accelerated  debris 
ions.  Finally,  Figs.  2c  and  2d  show  the  ions  and  fields  at  t  = 
0.0617  and  t  =  0.0771  respectively.  The  general  features  are 

similar  to  t  =  0.0411  (Fig.  2b).  There  are  intense  turbulent  fields  at  the 
debris-air  Interface  which  transfer  momentum  from  the  debris  to  air  ions, 
and  also  thermalize  the  particles.  Debris  ions  extend  out  to  x  >  10.0  rj^^ 
by  t  =  0.0617  as  well  as  the  B  field  and  low  level  turbulence.  In 

this  run  there  is  no  indication  of  Larmor  coupling.  One  interesting 
feature  seen  in  Fig.  2d  is  the  "structuring"  of  the  magnetic  pulse.  In 
particular,  there  is  a  sharp  gradient  that  has  developed  at  x  ~  4.8  r-j^^ 
which  produces  a  laminar  electric  field.  This  field  acts  to  energize  air 
ions  to  such  that  V  ~  1.8  at  x  ~  5.0  r^^. 

IV.  DISCUSSION 

We  have  presented  the  results  of  simulations  of  very  early  time 
processes  (l.e.,  t  <  !7^^)  using  a  ID  hybrid  code.  The  purpose  of  this  work 
has  been  to  elucidate  important  physical  phenomena  rather  than  to  perform 
an  actual  early  time  simulation  of  a  HANF.  The  focus  of  this  study  has 
been  primarily  on  coupling  mechanisms;  an  area  that  has  been  somewhat 
controversial  within  the  DNA  community.  We  have  found  that  there  is  no 
"single"  coupling  mechanism  but  that  several  mechanisms  can  be  effective  at 
various  times  throughout  the  debris  expansion  phase.  In  the  first 
simulation,  initialized  with  a  "broad"  magnetic  pulse,  the  initial  coupling 


occurs  through  laminar  electric  field  at  the  leading  edge  of  the  pulse. 
Subsequent  to  this,  turbulent  coupling  occurs  at  the  initial  debris~air 
interface  while  Larmor  coupling  occurs  within  the  main  body  of  the  magnetic 
pulse.  In  the  second  simulation,  initialized  with  a  "narrow"  magnetic 
pulse,  the  only  coupling  observed  was  because  of  the  electric  field 
(laminar  and  turbulent). 

A  very  important  result  of  this  study  is  the  generation  of  very 
energetic  debris  ions.  This  is  caused  by  the  laminar  electric  field  at  the 
leading  edge  of  the  magnetic  pulse.  A  fraction  of  the  debris  ions  can  be 
accelerated  up  to  roughly  three  times  the  initial  debris  expansion  velocity 
(i.e.,  Vji  ~  well-known  that  the  initial  formation  of  ion 
debris  patches  occurs  on  a  time  scale  faster  than  that  calculated  from  the 
initial  expansion  velocity.  Thus,  the  accelerated  debris  ions  observed  in 
the  simulations  is  consistent  with  HANE  observations. 

Finally,  this  study  dramatically  highlights  the  importance  of  kinetic 
phenomena  for  the  understanding  of  early  time  processes.  Kinetic  effects 
such  as  plasma  turbulence,  reflected  and  accelerated  ions,  and 
therraalization  of  the  ions  are  not  fluid  processes,  and  cannot  be  recovered 
using  only  MHD  codes.  Detailed  studies  are  necessary  using  codes  such  as 
the  one  described  in  this  paper  to  properly  understand  the  dynamics  of 
early  time  HANE  expansions. 
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Fig.  1  —  Evolution  of  the  ambient  magnetic  field  (bottom  curve),  debris  and  air  ions  (middle 
curves/dots),  and  electric  field  in  x-direction  (top  curve)  as  a  function  of  time.  The  parameters  are 
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Fig.  1  (Cont’d)  —  Evolution  of  the  ambient  magnetic  field  (bottom  curve),  debris  and  air  ions 
(middle  curves/ dots),  and  electric  field  in  x-direction  (top  curve)  as  a  function  of  time.  The 
parameters  are  =  4,  /.g  —  Iru,  and  njtta  —  16.  (a)  t  =  0.0059  (b)  t  =  0.0386 

nr';  (c)  t  =  0.0874  nr';  and  (d)  t  =  0.1342  nr'. 
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Fig.  1  (Cont’d)  —  Evolution  of  the  ambient  magnetic  field  (boi.om  curve),  debris  and  air  ions 
(middle  curves/dots),  and  electric  field  in  x-direciion  (top  curve)  as  a  function  of  lime.  The 
parameters  are  =  VJV^  -  4,  —  2ri,,,  and  nj “  16.  (a)  t  =  0.0059  fi,"';  (b)  t  -  0.0386 

n,-';  (c)  t  =  0.0874  n,-';  and  (d)  t  *  0.1342  n,"'- 
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Fig.  1  (Cont’d)  —  Evolution  of  the  ambient  magnetic  field  (bottom  curve),  debris  and  air  ions 
(middle  curves/dots),  and  electric  field  in  x-direction  (top  curve)  as  a  function  of  time.  The 
parameters  are  «»  VJV^  =  4,  —  2ri,,  and  nj n„  -  16.  (a)  t  =  0.0059  fl”';  (b)  t  *  0.0386 

nr'-,  (c)  t  =  0.0874  nr';  and  (d)  t  =  0.1342  n,"'. 
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Fig.  2  (Cont'd)  —  Evolution  of  quantities  described  in  Figure  1  but  with  —  VJ  = 
Lb  ~  0.5  ru,  and  nj -  64.  (a)  t  =  0.0206  n,"';  (b)  t  =  0.0411  n,-’:  (c)  t  =  0.0617  0," 
and  (d)  t  =  0.0771  n,"*. 
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APPENDIX 

HYBRID  CODE  DESCRIPTION 


We  describe  here  the  numerical  details  of  the  hybrid  code  used  in  this 
paper.  This  description  is  taken  from  Winske  and  Leroy  (1983)  with  minor 
changes  necessary  due  to  the  HANE  application.  Some  of  the  techniques 
(e.g.,  particle-in-cell  method)  are  standard  and  well  known,  while  others 
(e.g.,  solution  of  the  field  equations)  have  been  applied  primarily  to 
fusion-related  problems  and  hence  are  generally  not  familiar  to  the  space 
science  community.  The  discussion  is  divided  into  five  parts:  (A)  an 
overall  description  of  the  simulation  methods,  (B)  treatment  of  the  ions, 
(C)  solutions  of  the  field  and  electron  fluid  equations,  (D)  refinements  to 
the  basic  model,  and  (E)  discussion  of  the  initial  and  boundary 
conditions.  The  methods  pertaining  to  the  first  four  parts  are  rather 
general  and  could  thus  be  applied  to  any  problem  where  a  'hybrid'  (kinetic 
ion,  fluid  electron)  description  is  appropriate.  Specific  adaptation  of 
the  method  to  the  HANE  study  enters  primarily  through  the  last  topic,  the 
initial  and  boundary  conditions. 

A.  Overall  Description 

The  simulation  method  is  particle-in-celi  (Morse,  1970)  with  the 
electrons  treated  as  a  fluid  (l.e.,  cell)  quantity  rather  than  as  discrete 
particles.  The  assumption  of  a  massless,  charge  neutral  electron  fluid 
thus  eliminates  the  restriction  of  time  and  spatial  lengths  to  he  inverse 
electron  plasma  frequency  and  electron  Debye  length,  usually  associated 
with  full  particle  codes.  The  hybrid  models  discussed  here  were  originally 
developed  for  and  extensively  applied  to  dynamical  studies  of  high  density 
pinch  experiments  by  Chodura  (1975),  Sgro  and  Nielson  (1976)  and  Hamasakl 
et  al.  (1977).  These  models  apply  strictly  to  the  case  of  only  one  spatial 
dimension;  extension  to  two  dimensions  is  possible,  although  nontrivial,  as 
discussed  later.  Even  with  only  one  spatial  dimension  (which  is  along  the 
x-axis  and  thus  implies  3/ 3y  =  3/3z  =  0),  all  three  velocity,  magnetic 
field  and  electric  field  components  are  included.  The  simulation  region 
has  length  L,  divided  into  N  cells,  each  of  length  Ax  =  L/N.  There  is  an 
additional  ('ghost')  cell  at  each  end  of  the  system,  which  helps  in  keeping 
track  of  the  particles  entering  and  exiting  the  system  and  in  setting  up 


the  boundary  conditions.  The  ceil  quantities  (i.e.,  fields,  electron  fluid 
properties,  ion  velocity  moments)  are  specified  at  the  cell  centers. 
Because  the  restriction  to  electron  spatial  and  temporal  scales  has  been 
eliminated,  much  larger  time  steps  can  be  used.  Typically,  the  time 
step  (At)  is  limited  by  the  condition  that  ton  gyromotlon  is  well  defined 
and  the  cell  size  is  constrained  by  the  condition  that  the  fastest  ions  do 
not  traverse  one  ceil  in  one  time  step.  The  cell  size  and  system  length 
are  also  chosen  to  resolve  length  scales  of  interest  to  the  problem. 

3.  Dynamics 

The  ion  component  is  modeled  by  a  discrete  set  of  particles.  The  ion 
distribution  is  advanced  in  time  by  stepping  forward  each  particle  in  time 
under  the  Influence  of  the  local,  self-consistent  Lorentz  force.  The 
motion  of  the  ions  in  a  four-dlmens tonal  phase  space  (v^,  v^,  v^,  x)  is 
solved  by  the  particle-ln-cell  technique,  using  a  second-order-accurate  but 
non-reversible  scheme  (Nielson  and  Lewis,  1976).  The  equations  for  the  ion 
advance  (with  superscripts  denoting  the  time  level)  are: 


v"^^^  +  (h/2) 


+  h(E^  +  + 


.  At  vl/2 


— 0  _ 

V  X  B 


+ 


where  h  =  At  e/ra^,  f  =  l-(h^/2)B*^  •  B^,  g  =  rh/2)(v  •  B^),  e  is  the 

charge  and  m^  is  the  mass  of  the  tons.  The  electric  [e)  and  magnetic 

fields  (b)  are  evaluated  at  the  particle  position  x^.  P  is  a  mean  friction 

force,  P  =  -eg  •  J,  exerted  by  the  electrons  as  a  macroscopic  force 

only,  J  is  the  current  and  g  represents  a  phenomenological  anomalous 

resistivity  which  gives  rise  to  Ohmic  heating,  as  will  be  explained 

later.  After  all  of  the  ions  have  been  advanced,  the  ion  density  fn^^)  is 

updated,  by  averaging  over  the  positions  of  ail  the  particles  in  each 

cell.  The  ton  mean  velocities  fv .  1  are  also  needed  at  time  level  1.  The 

^~1 ' 

velocities  are  pushed  ahead  one  half  time  step, 

-1  1/2  1  1 

v^  =■  v^'^  +  (■h/2)fE^  +  - +  P^),  (A2) 


is  the 
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using  the  fields  at  the  new  particle  positions  x  ,  in  order  to  accumulate 
the  velocity  moments  V^. 

The  assumption  of  charge  neutrality  then  implies 

n  =  n  =  n,  (A3) 

e  i 


and  the  continuity  equation  gives 


3t 


(n  +  n.  )  = 
e  i 


I"  J 

3x  X 


~  enfv  -  V  ) 
3x  ix  eX' 


(AA) 


or 


V  =  V,  . 
ex  ix 


(A5) 


Thus,  the  electron  density  and  one  component  of  the  electron  fluid  velocity 
are  determined. 


C.  Field-Fluid  Equations 


With  the  assumption  of  zero  electron  mass  the  electron  momentum 
equation  reduces  to 


dV  V  X  B  -  . 

n  ra  =  0  =  -en  (E  +  — - )  -  -r—  p  x  -  n  P 

r  ^  3x  e 


e  e  dt 


e~ 


(A6) 


where  p^  is  the  scalar  pressure  and  is  the  velocity  of  the  electron 
fluid.  The  transverse  (y  and  z)  components  of  (A6)  can  be  written  as  Ohm's 
law 


J  =  2”^  * 


V  X  B 

::s — 2), 


(A7) 


In  general,  the  magnetic  field  is  inclined  at  some  angle  with  respect  to 

the  X  axis;  3  =  iBlTb  ,  b  ,  b  1  with  b^  +  +  b^  =  1,  and  if  a-fn.) 

~  x  y  z'^  x  V  z  Hi 


X  y  z "  x  y  z 

represents  the  conductivity  parallel  (perpendicular)  to  the  magnetic  field. 


the  conductivity  tensor  is  (Krall  and  Trlvelpiece,  1973): 


I 


a,  a  a 

1  y  z 

a  0-  a  , 

y  2  X 

o  a  a- 

z  X  3 


where 


0,  +  f0  -  0, lb  b 

i  (I  1-  X  X 


0.  +  f0-0,lbb 

1  I!  i  -  y  y 


0,  +  (0.  -  0,1b  b 

i  ^  II  1''  z  z 


0  =•  (0,,  -  o,  lb  b 

X  '■II  1-^  y  z 

0  =  (0,  -  0 ,  lb  b 

y  '■II  i-'  X  z 

0  =  (a ,  -  0,lbb. 

z  II  I-*  X  y 

The  transverse  fields  can  be  expressed  in  terms  of  the  vector 
potentials  in  the  usual  manner  fs  “  £  A;  E  =  -  fl /c )  f 3A/ 9t  1 1 : 

9A 

%  ■  ^ 


z  3x 


,  3A 
1  z 


and  7  •  3  =  0  implies 


B  =  constant. 

X 


(AlO) 


(The  last  field  component,  is  derived  from  the  x-component  of  (A6)  as 

will  be  shown  Later.) 


.  J 
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Neglecting  the  displacement  current.  Ampere's  law  can  be  expressed 

4ir 


-7 


3x 


J 

C  V 


9 

3“A 


4tT 


3x‘ 


- ^  J  . 

c  z 


(All) 


Substituting  in  Ohm's  law  (A7)  and  the  definitions  of  the  field 
components,  Eq.  (11)  becomes  a  coupled  system  of  equations: 
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3A 

3A 

3A 

3A 

— z  =  f 

^  3t  ’ 

z 

—  1 

2  1 
dX 

3t  ’ 

3x  ’ 

3x 

9 

3“A 

3A 

3A 

3A 

3A 

^  =  f 

'  3t  ’ 

z 

— 1 
3x  - 

,  2  ^2 
3x 

3t  ’ 

3x  ’ 

(A12) 


In  order  to  solve  this  system  a  fully  Implicit,  space  centered 


difference  scheme  is  used  (Richtmyer  and  Morton,  1967).  If  we  let  u^  =  Ay 

at  the  jth  ceil  position  and  n-th  time  level  (and  w^  =  A  in  the  same  way), 

I  2 


the  left  hand  side  of  (A12)  is  differenced  as  follows: 


3^  A 


n+l  ,  n+1  »  n+I 

u...  +  u.  ,  -  2u. 

.1+1  .1-t  ,1 


(A13) 


3x  Ax 

while  on  the  right  hand  side 
3A 


n+l  n 

u .  -  u. 

y  _  j _ j 


3t 


At 


and 


n+l  n+l 

1+1  1-1 

2  Ax 


(A14) 


'.Vhen  this  is  substituted  into  (A12),  the  resulting  system  of  equations  can 
be  written  symbolically  as 


n+1 

u. 

X  =  ■^ 

~j  n+l 
w. 

J 


(A16) 


and  Aj,  and  Cj  are  2x2  matrices  that  depend  only  on  components  of  a  and 
V  (evaluated  at  cell  i  at  time  level  n)  and  constants,  while  D  also 
depends  on  u?,  w^.  Thus,  all  the  coefficients  A,  B,  C  and  D  are  explicitly 
known.  Mote  that  there  are  N  cells  in  the  computation  mesh  j  =  2,  ...  M 
with  ghost  cells  (.1*0  and  1=«N+1)  at  each  end. 

The  set  of  equations  (A15)  is  easily  solved.  Assuming  a  solution  of 
the  form 


X.  =  E.  •  X._^,  +  F. 


1=0,  1,  ...  M 


(A17) 


and  substituting  it  into  (A15),  it  then  follows  that 


E.  =  -  [b.  +  c,  •  E.  ,  A. 

F.  =  [b.  +  C,  •  E.  ,  •  fo  -  C  •  F.,  ,  1 


(A18) 


Because  all  the  A.,  B.,  C,  and  D.  are  known,  if  E»  and  are  known  (from 
=3  =3  =1  ~3  “0  ~0 

boundary  conditions,  as  explained  later),  (A18)  can  be  solved  in  ascending 

order  to  obtain  E,,  F,  ;  E-,  F„,  ...,  E„,  F,,.  Then  if  X  can  be 

"I  ~1  =2  ~2  ’  “N  ~N  ~n+l 

determined  (again  from  boundary  conditions),  (A17)  can  be  solved  in 
descending  order  to  obtain  2^-1’  *’*'  5o’ 

Once  the  components  of  the  vector  potential  are  known,  the  new  values 
of  the  electric  and  magnetic  field  follow  from  (a9)  and  the  components  of 
the  current  from  (All),  again  using  the  differencing  in  (A13)  and  (A14). 
The  transverse  components  of  V  are  then  easliv  found  from  the  current  and 
the  Ion  velocity  moments: 


V  .  =  V,  -  J./en 
e3  11  3 


1  =  y  or  z. 


(A19) 


The  last  electron  fluid  quantity  to  be  calculated  is  the  temperature. 
The  electron  energy  equation  can  be  written  as 


t  -k  ("Je)  ♦  h  f|  *  Ve  -a-f  ■  0- 


(A20) 


The  source  term  Q  includes  resistive  heating  (nJ  )  and  loss  mechanisms, 
such  as  thermal  conduction  or  radiation,  depending  on  the  application. 
Letting 


P  "  n  T  , 
^e  e  a’ 


(A21) 


(A20)  can  be  rewritten  in  the  form  (with  y  »  5/3) 


^  “  “  YP 

'3t  ex  Sx-^^e  ^ 


e  3x 


+  (y-1)Q 


(A22) 


The  differencing  for  3p  /3t  and  3V  /3x  is  identical  to  (A14);  for 

e  ex 

stability  reasons  "donor  cell"  differencing  is  used  for  the  convective  term 
(Rlchtmyer  and  Morton,  1967): 


ex  3x 


/  n+1  n+1  ■> 

f".)  ■  ”.,1-2) 

/•  n+1  n+1  -v 

’■Pe.l+l  ■  Pe.1  ^ 


V  .  >  0 
exj 


V  ^  <  0 
exj 


Again,  a  tridiagonal  system  is  obtained. 


(A23) 


(A24) 


except  that  now  all  quantities  are  simple  scalars.  As  before,  an  assumed 
solution  of  the  form 


(A25) 


leads  to 


‘J  -  -  [»3  -  CjEj.,  r‘A^, 


Appropriate  boundary  conditions  lead  to  a  determination  of  Eq,  Fq  (and  then 
to  all  Ej,  Fj  in  ascending  order,  using  (A26))  and  and  thus  to  all  Xj 
in  descending  order  through  (A25)).  Since  n^^)  is  known,  is  then 
obtained  from  (A21).  An  alternative  to  (A22)  which  is  sometimes  useful 
(Sgro,  1978)  is  to  write  the  energy  equation  (A20)  in  terms  of  the 
entropy  [s  “  ln(p^n^^)1  Instead  of  the  pressure. 

Finally,  the  x  component  of  the  electron  momentum  equation  (A6)  can  be 
solved  for 


"X" 


-  1  (V  X  B  )  -  i-  ^  . 
c  ^~e  ne  3x 


(A27) 


This  electric  field  is  needed  to  maintain  charge  neutrality. 

The  computation  loop  through  one  time  step  can  thus  be  summarized  as 
follows; 

1.  Advance  the  ions  one  time  step  and  then  calculate  the  ion  velocity 

moments  (n^  =•  n^,  =>  V^y, 

2.  Compute  the  plasma  conductivity  g,  according  to  some  prescription. 
(The  resistivity  g  is  g  ). 

3.  Solve  the  coupled  equations  for  A„,  A_;  we  can  then  easily 

y  2 

calculate  Ey,  E^,  By,  B^  (B^^  “  constant).  From  the  currents 
(Ampere's  law)  and  the  ion  moments  calculate  the  other  two 
electron  velocity  components,  V^y  and 

4.  The  electron  temperature  is  obtained  next  from  the  solution  of  the 
differenced  energy  equation. 

5.  is  then  calculated  from  the  x-component  of  the  electron 
momentum  eauatlon. 

6.  Thus,  all  field  and  electron  fluid  components  are  known  and  we  are 
ready  to  move  the  ions  again. 


D.  Refinements 


The  simulation  scheme  described  thus  far  Is  rather  general,  the  only 
assumptions  being  quasineutrality,  zero  electron  mass  and  one  spatial 
dimension.  All  of  these  conditions  can  be  replaced  by  more  appropriate  j 

ones,  as  the  physical  situation  dictates.  For  example,  for  the  study  of 
low  frequency  ion  waves  an  adiabatic  electron  model  Is  more  appropriate; 
such  a  model  has  been  successfully  used  by  Okuda  et  al.  (1978).  In  some 
situations  phenomena  in  the  lower  hybrid  frequency  range  are  of  Interest. 

In  this  case  electron  Inertia  effects  are  non-negligible;  they  have  been 
included  in  the  one-dimensional  model  of  Liewer  (A1976)  and  the  two- 
dimensional  model  of  Hewett  and  Nielson  (A1978).  Two-dimensional 
simulation  models  with  m^  =  0  have  been  successfully  employed  by  Byers  et 
al.  (1978),  Hewett  (1980)  and  Harned  (1982). 

The  HANE  simulations  described  in  this  paper  take  the  conductivity  to 
be  constant.  A  further  level  of  sophistication  is  to  include  anomalous 
processes  due  to  the  mlcrophyslcs,  which  is  occurring  on  time  and  distance 
scales  shorter  than  those  resolved  in  the  hybrid  model,  by  means  of  more 
complicated  transport  coefficients.  This  is  done  by  expressing  the 
conductivity  (or  resistivity)  as  a  sum  of  two  terms,  one  representing 
classical  (Spltzer)  effects  and  the  other  due  to  anomalous  effects  arising 
from  microinstabilities  due  to  cross-field  currents.  Two  types  of 
anomalous  transport  coefficients  have  been  used  successfully  to  model  the 
behavior  of  laboratory  plasmas:  one  type  uses  a  seml-emplrical  expression 
(Chodura,  1975;  Sgro  and  Nielson,  1976),  the  other  is  based  on  a 
quasllinear  analysis  of  known  instabilities  (Davidson  and  Krall,  1977). 

The  transport  processes  can  be  further  refined  to  include  electron  thermal 
conductivity  (Sgro,  1980)  and  anomalous  ion  heating  (Hamasaki  et  al., 

1977).  Multispecies  ions  can  also  be  included  in  a  straightforward  manner 
(Sgro,  1980;  Sgro  and  Winske,  1981). 


E.  Initial  and  Boundary  Conditions 

In  order  to  simulate  the  piston  formation  stage  of  the  HANE  problem, 
appropriate  initial  and  boundary  conditions  have  been  implemented  in  the 
hybrid  code.  We  assume  initially  that,  a  dense  cool  debris  plasma  flows 
across  the  magnetic  field  into  a  tenuous  stationary  air  plasma.  In  this 
paper,  we  assume  the  debris  and  air  to  consist  of  the  same  ion  species. 
This  assumption  is  not  necessary  and  we  are  presently  Investigating  the 
effects  of  different  debris  composition.  We  assume  further  that  there  is  a 
magnetic  compression  of  finite  width  imbedded  in  the  leading  edge  of  the 
debris  stream,  formed  as  a  result  of  the  preceding  magnetic  compression 
stage  of  the  coupling  la  simply  maintained  at  the  right  boundary.  Ions  of 
all  species  are  effectively  absorbed  at  the  boundaries.  The  final  boundary 
condition  is  on  the  magnetic  field  which  is  held  constant  at  the  ambient 
(air)  value  at  both  boundaries.  The  magnetic  compression  is  taken  to  be 
proportional  to  the  electron  compression.  The  width  of  the  magnetic 
compression  is  an  input  parameter  in  the  code.  Thus  initially  in  the 
simulation  space  we  establish  a  debris  ion  stream  extending  a  finite 
distance  from  the  left  boundary  into  the  uniformly  distributed  air  ions.  A 
step  function  compression  in  the  magnetic  field,  narrower  than  the  debris 
stream,  is  Initialized  at  the  leading  edge  of  the  debris.  The  magnetic 
field  is  initially  uniform  on  both  sides  of  the  compression.  During  the 
simulation  runs,  debris  ions  are  continuously  injected  from  the  left 
boundary.  This  is  accomplished  by  re-inltlalizlng  the  debris  ion 
distribution  in  the  left  ghost  cell  each  time  step. 


» 


I 


27 


REFERENCES 


Byers,  J.A.,  B.I.  Cohen,  W.C.  Condit ,  and  J.D.  Hanson,  "Hybrid  Simulations 
of  Quasineutral  Phenomena  in  Magnetized  Plasma,"  J  Comp.  Phys . ,  27, 
363,  1978. 

Chodura,  R.,  "Hybrid  Fluid-Partiole  Model  of  Ion  Heating  in  High-Maon- 
Nuraber  Shook  Waves,"  Nuol.  Fusion,  15,  55,  1975. 

Davidson,  R.C.,  and  N.A.  Krall ,  "Anomalous  Transport  in  High-Temperature 
Plasmas  with  Applications  to  Solenoidal  Fusion  Systems,"  Nuol,  Fusion, 
U*  131  3,  1  977. 

Hamasaki ,  S. ,  N.A.  Krall,  C.E.  Wagner,  and  R.N.  Byrne,  "Effect  of 
Turbulence  on  Theta  Pinch  Modeling  by  Hybrid  Numerical  Models,"  Phys . 
Fluids,  65,  1977 

Harned,  D.S.,  "Quasineutral  Hybrid  Simulation  of  Macroscopic  Plasma 
Phenomena,"  J.  Comp.  Phys.,  47,  452,  1982. 

Hewett,  D.W.,  "A  global  method  for  Solving  the  Electron-field  Equations  in 
a  Zero-inertia-electron-hybrid  Plasma  Simulation  Code,"  J.  Comp. 
Phys . ,  378,  1980. 

Hewett,  D.W.  and  C  .W .  Nielson,  "A  Multidimensional  Quasi  neutral  Plasma 
Simulation  Model,"  J.  Comp.  Phys.,  29,  219,  1978. 

Krall,  N.A.  and  A.W.  Trivelpiece,  Principles  of  Plasma  Physics,  p.  119, 
McGraw-Hill,  New  York  1973. 

Liewer,  P.C.,  "Numerical  Studies  of  Ion  Reflection  in  Collisionless  Theta- 
Pinch  Implosions  Using  a  Hybrid  Vlasov-Fluid  Model,"  Nucl.  Fusion,  1 6 , 
817,  1976. 

Longmire,  C.L.,  "Notes  on  Debris-Air-Magnetic  Interaction,"  The  Rand 
Corporation  Report  RM-3386-PR,  1963.  AD296597. 

Morse,  R.L.,  "Multidimensional  Plasma  Simulation  by  the  Particle-in-Cell 
Method,"  Methods  in  Computational  Physics,  Vol,9,  Ed.  by  B.  Alder,  S. 
Fernbach  and  M.  Rotenberg,  p.  213,  Academic  Press,  NY,  1970. 

Nielson,  C  .W .  and  H.R.  Lewis,  "Particle- code  Models  in  the  Nonradiatlve 
Limit,"  in  Methods  in  Computational  Physics,  Vol.  16,  edited  by  J. 
Killeen,  B.  Alder,  S.  Fernbach  and  M.  Rotenberg,  p.  367,  Academic 


Okuda ,  H.,  J.M.  Dawson,  A.T.  Lin  and  C.C.  Lin,  "Quasi-neutral  Parti  ole 
Simulation  with  Application  to  Ion  Wave  Propagation,"  Phys.  Fluids, 
476,  1978. 

Papadopoulos ,  X.,  R.C.  Davidson,  J.M.  Dawson,  I.  Haber,  D.A.  Hammer,  N.A. 
Krall ,  and  R.  Shanny,  "Heating  of  Counterstreaming  Ion  Beams  in  an 
External  Magnetic  Field,"  Phys.  Fluids,  14,  849,  1971. 

Ripin,  3.H.,  J.  Grun,  S.  Kaoenjar,  E.A.  McLean,  and  J.A.  Stamper, 
"Introduction  to  the  Laser-HANE  Experiment  and  Summary  of  Low-Pressure 
Interaction  Results,"  NRL  Memo  Report  5268,  1984.  ADI  38945. 

Sgro,  A.G.  and  C.W, Nielson,  "Hybrid  Model  Studies  of  Ion  Dynamics  and 
Magnetic  Field  Diffusion  During  Pinch  Implosions,"  Phys.  Fluids,  1 9, 
126,  1976. 

Sgro,  A.G.,  "Calculation  of  the  Effects  of  Incomplete  Preionization  in  High 
Voltage  Theta  Pinches,"  Phys.  Fluids,  21 ,  1410,  1978. 

Sgro,  A.G.,  "Simulation  of  the  ZT-S  reversed  Field  Pinch,"  Phys.  Fluids, 
1055,  1980. 

Sgro,  A.G.  and  D.  Winski ,  "Simulation  of  the  Formation  of  the  PS-1 
Spheromak,"  Phys.  Fluids,  24,  1156,  1981. 

Winske,  D.  and  M.M.  Leroy,  "Hybrid  Simulation  Tech  Applied  to  Earth's  Bow 
Shock,"  Comp.  Simulations  of  Space  Plasmas  -  Selected  Lectures  at  1st 
ISS,  D.  Reidel,  Holland,  1984. 


29 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 


ASSISTANT  SECRETARY  OF  DEFENSE 
COMM,  CMD,  CONT  7  INTELL 
WASHINGTON,  D.C.  20301 

DIRECTOR 

COMMAND  CONTROL  TECHNICAL  CENTER 
PENTAGON  RM  BE  685 
WASHINGTON,  D.C.  20301 
OICY  ATTN  C-6S0 
01CY  ATTN  C-312  R.  MASON 

DIRECTOR 

DEFENSE  ADVANCED  RSCH  PROJ  AGENCY 
ARCHITECT  BUILDING 
UOO  WILSON  BLVD. 

ARLINGTON,  VA.  22209 
OICY  ATTN  NUCLEAR 

MONITORING  RESEARCH 
OICY  ATTN  STRATEGIC  TECH  OFFICE 

DEFENSE  COMMUNICATION  ENGINEER  CENTER 
1860  WIEHLE  AVENUE 
RESTON,  VA.  22090 
OICY  ATTN  CODE  R410 
OICY  ATTN  CODE  R812 

DEFENSE  technical  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA .  2231A 

02CY 

DIRECTOR 

defense  NUCLEAR  AGENCY 
WASHINGTON,  D.C.  20305 
OICY  ATTN  STVL 
04CY  ATTN  TITL 
OICY  ATTN  DDST 
03CY  ATTN  RAAE 

COMMANDER 
FIELD  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
KIRTLAND,  AFB,  NM  87115 
OICY  ATTN  FCPR 


DEFENSE  NUCLEAR  AGENCY 
S AO/ONA 

BUILDING  20676 
KIRTLAND  AF9,  NM  87115 
OICY  D.C.  THORNBURG 

DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AFB,  NM  87115 

OICY  ATTN  DOCUMENT  CONTROL 

JOINT  CHIEFS  OF  STAFF 
WASHINGTON,  D.C.  20301 

OICY  ATTN  J-3  WWMCCS  EVALUATION 
OFFICE 

DIRECTOR 

JOINT  STRAT  TGT  PLANNING  STAFF 
OFFUTT  AFB 
OMAHA,  NB  68113 

OICY  ATTN  JSTPS/JLKS 
OICY  ATTN  JPST  G.  GOETZ 

CHIEF 

LIVERMORE  DIVISION  FLD  COMMAND  DNA 
DEPARTMENT  OF  DEFENSE 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 
OICY  ATTN  FCPRL 

COMMANDANT 

NATO  SCHOOL  (SHAPE) 

APO  NEW  YORK  09172 

OICY  ATTN  U.S,  DOCUMENTS  OFFICER 

UNDER  SECY  OF  DEF  FOR  RSCH  R  ENGRG 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  D.C.  20301 

OICY  ATTN  STRATEGIC  R  SPACE 
SYSTEMS  (OS) 

WWMCCS  SYSTEM  ENGINEERING  ORG 
WASHINGTON,  D.C.  20305 
OICY  ATTN  R.  CRAWFORD 


COMMANDER /DIRECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OUY  ATTN  DELAS-EO,  F.  NILES 

DIRECTOR 

BHD  ADVANCED  TECH  CTR 
HUNTSVILLE  OFFICE 
P.O.  BOX  1500 
HUNTSVILLE,  AL  35807 

01CY  ATTN  ATC-T  MELVIN  T.  CAPPS 
01CY  ATTN  ATC-0  W.  DAVIES 
Q1CY  ATTN  ATC-R  DON  RUSS 

PROGRAM  MANAGER 
BHD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

01CY  ATTN  DACS-3MT  J.  SHEA 

CHIEF  C-E-  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CMD 
PENTAGON  RM  1B269 
WASHINGTON,  D.C.  20310 

01CY  ATTN  C-  E-SERVICES  DIVISION 

COMMANDER 

FRADCOM  TECHNICAL  SUPPORT  ACTIVITY 
DEPARTMENT  OF  THE  ARMY 
FORT  MONMOUTH,  N.J.  07703 

OUY  ATTN  DRSEL-NL-RD  H.  BENNET 

01CY  ATTN  ORSEL-PL-EHV  H.  BOMKE 

OUY  ATTN  J.E.  QUIGLEY 

COMMANDER 

U.S.  ARMY  COHM-ELEC  ENGRG  INSTAL  AGY 
FT.  HUACHUCA,  AZ  85613 

OUY  ATTN  CCC-EHEO  GEORGE  LANE 


COMMANDER 

U.S.  army  NUCLEAR  AND  CHEMICAL  AGENCY 

7500  SAC<LIC<  ROAD 

BLDG  2073 

SPRINGFIELD,  VA  22150 
OUY  ATTN  library 

DIRECTOR 

U.S.  army  PALLISTIC  research 

LABORATORY 

ABERDEEN  PROVING  GROUND,  MD  21005 
01CY  ATTN  TECH  LIBRARY, 

EDWARD  BAICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 

FT.  MONMOUTH,  NJ  07703 

OUY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 

REDSTONE  ARSENAL,  AL  35809 
OUY  ATTN  JIM  gamble 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS 

ACTIVITY 

WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OUY  ATTN  ATAA-SA 
OUY  ATTN  TCC/F.  PAYAN  JR. 

OUY  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  electronic  SYSTEMS  COMMAND 

WASHINGTON,  D.C.  20360 

OUY  ATTN  NAVALEX  034  T.  HUGHES 
OUY  ATTN  PME  11^ 

OUY  ATTN  PME  117-T 
OUY  ATTN  CODE  5011 


COMMANDER  COMMANDING  OFFICER 

U.S.  ARMY  FOREIGN  SCIENCE  6  TECH  CTR  NAVAL  INTELLIGENCE  SUPPORT  CTP 
220  7TH  STREET,  NE  4301  SUITLAND  ROAD,  BLDG.  5 

CHARLOTTESVILLE,  VA  22901  WASHINGTON,  D.C.  20390 

OUY  ATTN  ORXST-SD  OUY  ATTN  MR.  DUBBIN  STIC  12 

OUY  ATTN  N  I  SC -50 

COMMANDER  OUT  ATTN  CODE  5404  J.  GALET 

U.S,  ARMY  MATERIAL  DEV  8  READINESS  CMD 

5001  EISENHOWER  AVENUE  COMMANDER 

ALEXANDRIA,  VA  22333  NAVAL  OCCEAN  SYSTEMS  CENTER 

OUY  ATTN  ORCLDC  J.A.  BENDER  SAN  DIEGO,  CA  97152 

OUY  ATTN  J.  FERGUSON 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  O.C.  20375 

OICY  ATTN  CODE  4700  S.  L.  Ossakou 
26  CYS  IF  UNCLASS.  1  CY 
IF  CLASS) 

01CY  ATTN  CODE  4701  I  Vitkovitsky 
01CY  ATTN  CODE  4780  J.  Huba  (50 

CYS  IF  UNCLASS,  1  CY  IF  CLASS) 
01CY  ATTN  CODE  7500 
OICY  ATTN  CODE  7550 
01CY  ATTN  CODE  7580 
OICY  ATTN  CODE  7551 
OICY  ATTN  CODE  7555 
OICY  ATTN  CODE  4730  E.  MCLEAN 
OICY  ATTN  CODE  4108 
OICY  ATTN  CODE  4730  8.  RIPIN 

20CY  ATTN  CODE  2628 

COr"ANDER 

NAVAL  SPACE  SURVEILLANCE  SYSTEM 
DAHLGREN,  VA  22448 

OICY  ATTN  CAPT  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MO  20910 
OICY  ATTN  CODE  F31 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  O.C.  20376 
OICY  ATTN  NSP-2141 

OICY  ATTN  NSSP-2722  FRED  WIMBERLY 
COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22448 

OICY  ATTN  CODE  OF-U  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 
OICY  ATTN  CODE  465 
OICY  ATTN  CODE  461 
OICY  ATTN  CODE  402 
OICY  ATTN  CODE  420 
OICY  ATTN  CODE  421 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND /XPO 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 
OICY  ATTN  XPDQO 
OICY  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  OPR  HAROLD  GARDNER 
OICY  ATTN  LKB 

KENNETH  S.W.  CHAMPION 
OICY  ATTN  OPR  ALVA  T.  STAIR 
OICY  ATTN  PHD  JURGEN  BUCHAU 
OICY  ATTN  0  JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA  ARTHUR  H.  GUENTHER 
OICY  ATTN  NTYCE  1LT.  6.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
OICY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  *  ACO 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
OICY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DI V I S ION / YSE A 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
OICY  ATTN  YSEA 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

OICY  ATTN  DC  HR.  LONG 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/OC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  HA  01731 

OICY  ATTN  OCKC  MAJ  J.C.  CLARK 


COMMANDER 

FOREIGN  TECHNOLOGr  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  A5433 
OlCr  ATTN  NICD  LIBRARY 

01CY  ATTN  ETOP  8,  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13441 

01CY  ATTN  DOC  LIBRARY/TSLO 
01CY  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01CY  ATTN  S2J 

STRATEGIC  AIR  COHMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01CY  ATTN  AOWATE  HAJ  BRUCE  BAUER 
01CY  ATTN  MRT 

01CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/ SK 
P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
H.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01 CY  ATTN  HNNL 

COMMANDER 

ROHE  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  O.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  S7115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 


EGSG,  INC- 

LOS  alamos  DIVISION 

P.O.  BOX  809 

LOS  alamos,  N''  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

01CY  ATTN  DOC  CON  FOR  U-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J  . 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

OICY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

01  CY 

ATTN 

0. 

SAPPENFIELD 

SANOIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

01CY  ATTN  DOC  CON  FOR  W.  BROWN 
OICY  ATTN  DOC  CON  FOR  A. 

THORNBROUGH 

OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

01CY  ATTN  DOC  CON  FOR  0.  DAHLGREN 

OICY  ATTN  DOC  COM  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  HURPHEY 

OICY  ATTN  OOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 


INSTITUTE 

FOR 

TELECOM  SCIENCES 

NATIONAL 

TELECOMBNICATIONS  S 

INFO 

ADMIN 

BOULDER, 

CO 

80303 

OICY 

ATTN 

A.  JEAN  (UNCLASS 

ONLY) 

OICY 

ATTN 

W.  UTLAUT 

01  CY 

ATTN 

D.  CROMBIE 

OICY 

ATTN 

L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01CY  ATTN  R.  GRUBB 
01CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
01CY  ATTN  I.  GARFUNKEL 
01CY  ATTN  T.  SALMI 
01CY  ATTN  V-  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  HOAD 
BURLINGTON,  MA  01303 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01CY  ATTN  L.  SLOAN 
01CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  9A701 

01CY  ATTN  J,  WORKMAN 
01CY  ATTN  C,  PRETTIE 
01CY  ATTN  S,  BRECHT 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  9812A 

01CY  ATTN  G,  KEISTER 
01CY  ATTN  D.  MURRAY 
01CY  ATTN  G,  HALL 
OlCt  ATTN  J.  KENNEY 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMFIRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
01CY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 

linthicum  road 
CLAR<saURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 

OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  technologies,  INC. 

606  WiLshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 
OICY  ATTN  R.  LELEVIER 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
OICY  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  ATTN  G.  MILLMAN 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OE  ALASKA 
FAIRBANKS,  AK  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 

OICY  ATTN  T.N,  DAVIS  (UNCLASS  ONLY) 
01 CY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVAN  I  A,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
OICY  ATTN  J.H.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARO 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  8  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

J  AYCOR 


KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MFAGHER 

KAMAN  TEMPO-CEN’-ER  FOR  ADVANCED 
STUDIES 

816  SYATE  STREET  (P.O  DRAWER  (3Q) 

SANTA  BARBARA,  CA  93107 
OICY  ATTN  OASIAC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  S.  GAMHILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO,,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  S  SPACE  CD,,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
'.u.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 


11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L,  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  HO  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEHRA 
OICY  ATTN  JOHN  OASSOULAS 


M.I.T.  LINCOLN  LABORATORY 
P.O,  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  L.  LOOGHLIN 
OICY  ATTN  0.  CLARK 


MCDONfJEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01 CY  ATTN  GEORGE  MROZ 
01CY  ATTN  U.  OLSON 
01CY  ATTN  R,W.  HALPRIN 
01CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


01  CY 

ATTN 

P.  FISCHER 

01  CY 

ATTN 

W.F.  CREVIER 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01  CY 

ATTN 

R.  BOGUSCH 

01  CY 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RALPH  KILB 

01  CY 

ATTN 

DAVE  SOWLE 

01  CY 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEI8E 

01  CY 

ATTN 

CONRAD  L.  LONGMIRE 

01  CY 

ATTN 

B.  WHITE 

01  CY 

ATTN 

R.  STAGAT 

MISSION 

RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUOUERQUE,  NEW  MEXICO  87106 
01CY  R.  STELLINGWERF 
01CY  H.  ALME 
01CY  L.  WRIGHT 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01CY  ATTN  IRVING  L.  KOESKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND ,  CA  94610 
ATTN  A.  THOMSON 

R  ft  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


01  CY 

ATTN 

FORREST  GILMORE 

01  CY 

ATTN 

WILLIAM  B.  WRIGHT 

01CY 

ATTN 

WILLIAM  J.  KARZAS 

01  CY 

ATTN 

H.  ORY 

01CY 

ATTN 

C.  MACDONALD 

01CY 

ATTN 

R.  TURCO 

01CY 

ATTN 

L .  D  e  R  A  N  D 

01CY 

ATTN 

W.  TSAI 

MI'RE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 


A  = 


37 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
C1CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOOO.RIDGE  DR. 

MCLEAN,  VA  22102 
01CY  J.  COCKAYNE 
01CY  E,  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
C1CY  ATTN  J.  CASPER 
01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNOOA 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  6.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
01CY  ATTN  RAY  L.  LEADA0RAND 
Q1CY  ATTN  6.  CARPENTER 
Q1CY  ATTN  G.  PRICE 
01CY  ATTN  R.  LIVINGSTON 
01CY  ATTN  V,  GONZALES 
01CY  ATTN  D,  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

01CY  ATTN  JOEL  GARBARINO 

TRW  DEFENSE  8  SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90273 


01  CY 

ATTN 

R. 

K.  PLEBUCH 

01  CY 

ATTN 

S. 

ALTSCHULER 

01CY 

ATTN 

D. 

DEE 

01  CY 

ATTN 

0/ 

STOCKWELL 

SNTF/1S75 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01303 
01CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01CY  ATTN;  N,  ZABUSKY 

DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 
02CY 


